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ABSTRACT 

The baryon acoustic oscillations (BAO) prior to recombination should be imprinted onto the 21cm 
emission background from the epoch of reionization through the underlying density perturbations. Us- 
ing an analytical approach for both matter power spectrum (CDM+baryons) and reionization process, 
we demonstrate the BAO induced signatures on the power spectrum of 21cm emission fluctuations. 
Future low-frequency radio telescopes such as LOFAR and MWA should be able to detect these weak 
BAO wiggles with an integration time of ~ 1 year. A combination of the BAO measurements at dif- 
ferent redshifts z w 1000 (CMB), z rs 10 (epoch of reionization) and z rs (clustering of galaxies) may 
allow one to set more robust constraints on the determinations of cosmological parameters including 
dark energy and its equation of state. 

Subject headings: cosmology: theory — large-scale structure of universe — diffuse radiation — inter- 
galactic medium 



1. INTRODUCTION 

Prior to recombination, free electrons coupled the 
baryons and photons tightly through Compton scatter- 
ing, and these three species moved together as a single 
fluid. In this relativistic plasma, the primordial small- 
scale perturbations propagated as sound waves, result- 
ing in the pressure-induce oscillations. The neutral gas 
can still retain some memory of such acoustic oscillations 
even after recombination, manifesting themselves in the 
last scattering surface seen as the harmonic series of max- 
ima and minima on the cosmic microwave background 
(CMB) at redshift z ss 1000. The longest wavelength 
of BAO with A « 100 Mpc imprinted on the large-scale 
structures is still visible in the local universe through 
the survey of 3D galax y distributions (jEisenstein et al. I 
PiOOfit [C^ole et al. Il2005h . 

Because BAO can be served as an ideal cosmic ruler for 
many cosmologic al applications espec i ally for the probe 
of dark energy llHu fc White! 11995 iBarkana fc LodTl 
120051 IWvithe. Loeb fc Geil Il2007f ). it is of great inter- 
est to explore how BAO evolve with cosmic time, and in 
particular how and when the peaks of BAO at smaller 
scales are washed out by emergence of larger structures. 
Indeed, in addition to the detections of BAO signatures 
in CMB at z ~ 1000 and galaxy spatial distributions 
at lower redshifts, one may be able to extract valuable 
information on BAO at redshifts around z « 10 from 
the study of 21 cm absorption/emission generated in the 
dark ages and epoch of reionization. This will comple- 
ment our knowledge of structure formation at this im- 
portant phase when galactic dark halos started to de- 
velop. The existing modes and positions of BAO at the 
epoch of reionization would make a sensitive diagnosis of 
nonlinear structures evolved by that time. Most impor- 
tantly, the statistical uncertainties in the determination 
of cosmological parameters including dark energy and its 
equation of state can be significantly reduced when more 
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independent measurements of BAO through the epoch of 
reionization are incorporated with the BAO features al- 
ready detected in CMB and large-scale structures of the 
local universe. Note that the 21 cm absorption/emission 
observations provide a tomographic imaging of the uni- 
verse at the epoch of reionization and dark ages, yielding 
many independent constraints on the theory of cosmol- 
ogy- 

The goal of this letter is to demonstrate the signa- 
tures of BAO on the redshifted 21 cm fluctuations and 
discuss the feasibility of detections with existing and 
planned low frequency telescopes such as 21CMA, LO- 
FAR and MWA. We focus on the 21cm emission of neu- 
tral hydrogen during the process of re ionization instead 
of th e 21 cm absorption at z > 20 (|Barkana fc Loebl 
2005). For the latter, the detection of BAO signatures 
turns to be much more difficult because of their longer 
wavelength and limitation of angular resolutions with 
exiting and even future radio telescopes. A sophisti- 
cated treatment of the problem requires the detailed 
knowledge of the history of reionization and radiative 
transfer of ionizing photons through gas density field, 
which may be achievable by numerical simulations (e.g. 
iFurlanetto. Sokasian fc Hern auist 20041 ). Here we would 
rather employ an analytic approach to simplify the prob- 
lem and highlight the essentials of physical process at the 
epoch of reionization. Throughout the paper we adopt 
a concordance cosmology of flo = 0.265, f^A = 0.735, 
Q b = 0.044, h = 0.71, n s = 1 and a 8 = 0. 772, as revealed 
by th e WMAP three-year observations (Sp ergel et al. I 
12001 . 

2. MATTER POWER SPECTRA 

We begin with the linear matter power spectrum 
P]i n (z,k) cx D 2 (z)k ns T 2 (k), which describes how the 
initial matter power spectrum fc™ 3 is modulated by the 
transfer function T(fc) and growth factor D(z). The 
baryon content, and thereby BAO, is incorporated in 
T(k) which can be approximately separated into the cold 
dark matter (CDM) and baryon components: T(k) = 
(n b /n )T b (k) + (n c /n )T c (k), where n c is the CDM den- 
sity parameter relative to the critical density at present, 
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Fig. 1. — Linear matter power spectra at z = 6 for three different 
matter contents: pure CDM (dot-dashed line), pure baryons (grey) 
and mixed CDM+baryons (solid line). The BAO features can be 
seen at k ~ 0.1 Mpc -1 in the mixed matter model. 



O c = Hq — fif,. We adopt the asymptotic solutions to 
both Th(fc) and TJk) n ear the sound horizon given by 
lEisenstein fe Hu I |l998), in which the suppression effect 
of baryons on scales below the sound horizon is included. 
Nonetheless, the linear matter power spectrum becomes 
inaccurate at smaller scales and later cosmic time. We 
employ a halo model to evaluate the nonlinear power 
spectrum (see ICoorav fc Sheth I I2002T) . which accounts 
for contributions from the single halo term Pih plus the 
clustering term P2h- We take the Press-Schechter formal- 
ism for the m ass function of dark halos and th e Navarro- 
Frenk- White (jNavarro. Frenk fc White1ll996| ) profile for 
the matter distribution inside each halo. Employment 
of the nonlinear matter power spectrum instead of the 
linear one is crucial in the sense that BAO can be erased 
entirely on scales of nonlinear structures, although this 
may not be a serious problem at the epoch of reioniza- 
tion. 

In Fig. 1 , we demonstrate the linear matter power spec- 
tra Pn n at redshift z = 6 for three matter contents: pure 
CDM, pure baryons and mixed CDM+baryons. Inclu- 
sion of baryons in the pure CDM model gives rise to 
the weak yet visible wiggles in the matter power spec- 
trum at wavenumber fc ~ 0.1 Mpc -1 . Next, we calculate 
the nonlinear matter power spectrum instead of Pi; n but 
leave the baryon content unchanged, and the result for 
z — 6 is shown in Fig. 2. The nonlinear matter power 
spectrum at large scales fc < 1 Mpc -1 remains roughly 
the same as Pn n , in which the weak BAO are clearly pre- 
sented. The prominent nonlinear structures dominate 
the matter power spectrum only at short wavelengths of 
fc > 1 Mpc -1 . It turns out that by the end of cosmic 
reionization at z w 6, BAO are still unaffected by the 
formation of nonlinear structures. This arises because 
the gravitationally bound systems such as dark halos and 
their associated large-scale structures at z — 6 have sizes 
much smaller than the typical scales (~ 100 Mpc) of 
BAO. In other words, many of the interesting modes of 
BAO should leave their imprints on the matter power 
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Fig. 2. — The matter power spectra at z = 6 predicted by linear 
density perturbation (dotted line) and halo model (solid line). For 
the latter, contributions of 1-halo and 2-halo terms are also dis- 
played. Linear theory breaks down at small scales with wavenum- 
ber beyond k > 1 Mpc" 1 , and BAO are unaffected by nonlinear 
structures by z = 6. 

spectrum before z — 6. 

3. 21CM POWER SPECTRA 

BAO signatures enter into the redshifted 21cm emis- 
sion background from the epoch of reionization through 
the underlying matter density fluctuations (<5). If we 
restrict ourselves to the 21cm emission generated from 
the neutral hydrogen in the surroundings of the ionized 
bubbles of first-generation luminous objects, the surface 
brightness of the emission can be evaluated through (cf., 
IZaldarriaga. Furlanetto fc Hernquist II2004J ) 



8T b t*T (l+5)(l-x) 



T = 16mKh- 1 




(1) 



where x = x e {l + 5 X ) is the ionization fraction, x e is the 
average ionization fraction and 5 X is the perturbation in 
the ionization fraction across t he sky, for which we will 
take t he reionization model of ISantos. Coorav fc Knoxl 
( 2005). The corresponding power spectrum of the 21cm 
emission can be written as 



P$t>(z, fc) = T 2 (1 - x e ) 2 Pss{z, fc) + x 2 e P Sm s m (z, k) 



2x e (l - x e )Pss 3: (z 1 k) 



(2) 



The three terms in the right-hand side represent the 
contributions of the matter power spectrum, the power 
from the perturbations in the ionization fraction, and the 
cross-correlation po wer, respective l y. For the latter two, 
we use the model of ISantos et al. I ([2003) to proceed our 
numerical computation. The angular power spectrum 
Ci{v) of the redshifted 21cm brightness sky is derived 
under the flat-sky approximation. 

The theoretically predicted 3D and 2D power spectra 
of the redshifted 21cm emission fluctuations are shown 
in Fig. 3 and Fig. 4, respectively. While the amplitudes 
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Fig. 3. — Power spectra of 21cm fluctuations at different rcdshifts 
from z = 30 to z = 6. The ionization fraction x e is also showed at 
each redshift. BAO wiggles occur at k ~ 0.1 Mpc -1 and nonlinear 
structures dominate at k > 1 Mpc -1 . The matter power spectrum 
at z = 6 is also plotted for comparison (dotted line). 




Fig. 4. — Angular power spectra of 21cm fluctuations at different 
redshifts from z = 30 to z = 6. BAO wiggles are presented at 
I ~ 500 — 3000, and their positions tend towards large i with the 
increase of redshift. 



of the 21cm power spectra themselves are relatively low, 
with a maximum value of ~ 10 mK, the BAO induced 
wiggles are clearly presented. It appears that an angular 
resolution of ~ 10 arcminute is needed in order to iden- 
tify the BAO features on the 2D power spectra, apart 
from the requirement of high sensitivity. Moreover, the 
positions of BAO tend towards large £ with the increase 
of redshift, indicating that radio arrays with baseline of 
~ 1 km will be needed to reveal these structures. We 
have also calculated the ID power spectrum of the 21cm 
emission along the line of sight but found that the BAO 
features are completely washed out due to project ion ef - 
fect, in agreement with the result of IWang fc Hu I (2006). 
Finally, it is pointed out that our analytical model does 



7r^ eg /129600 ac- 



not take the size distribution of ionized bubbles into ac- 
count. The typical sizes of the ionized bubbles near z rss 6 
can reach ~ 20 Mpc, which is already comparable to the 
BAO wiggles at large k or £. Whether or not the ionized 
bubbles produce oscillations on the same scales as BAO 
should be investigated in future study. 

4. DETECTABILITY 

Detection of the BAO signatures on the redshiftcd 
21cm fluctuations indeed poses a technique challenge for 
existing and planned low-frequency radio telescopes. We 
demonstrate the observability using 21 CentiMeter Array 
(21CMA, cosmo.bao.ac.cn), Low Frequency Radio Ar- 
ray (LOFAR, www.lofar.org) and Mileura Widefield Ar- 
ray (MWA, www.haystack.mit.edu/arrays/MWA), and 
only work with the angular power spectrum Ci at 
a fixed frequency. Supposing that strong fore- 
ground contamination can be entirely removed from 
the low frequency sky through cither the two- 
point correlation technique in frequency domain (e.g. 
IZaldarriaga. Furlanetto fc L lcrnaiiist 2003) or the pixel- 
by-pixel algorithm (jWang et al. I 1200(f ). we can esti- 
mate the variance in Ct through AC( = [2/(2.6 + 
1)/ S ky] 1/2 (C, + N t ), in which / sky = 

counts for the sky coverage, and N# = (u>/ s k y 
is the noise power spectrum if we adopt a Gaussian func- 
tion with width 9b for the experimental beam and use 
w^ 1 = 47T(7p ix /A r pi x to denote the contribution of the 
white noise with cr p i x and N v - lyi being the pixel noise 
and total number of pixels, respectively. In radio in- 
terferometric measurement, the pixel noise can be rep- 
resented in terms of brightness temperature as er p ; x = 

T S y S /r)v2NAvt, where T sys is the system temperature, r\ 
is the efficiency factor of telescope, N is the total number 
of independent baselines, Ais is the bandwidth, and t is 
the observing time. 

To proceed further, for 21CMA we take a system tem- 
perature of T sys = 250 K and an efficiency of i] = 0.64. 
The total dishes of 21 CM A are iVdish = 81 and the 
longest baseline is 6 km, which gives rise to an angu- 
lar resolution of 0b ~ 1 arcmin. We use a conservative 
value of 9i, = 2 arcmin in the present estimate. The sky 
coverage is, nevertheless, very small: / s k y — 10~ 3 . For 
LOFAR (compact core), the corresponding parameters 
are chosen to be: T sys = 100 K, rj = 0.64, iVdish = 32, 
9b = 3 arcmin, and / s k y ~ 0.1. We utilize the fol- 
lowing parameters for MWA: T sys = 200 K, rj = 0.64, 
^dish = 500, 9b — 5 arcmin, and / s k y ~ 0.4. The er- 
rors ACi in the measurement of the 21cm power spectra 
with 21CMA, LOFAR (core) and MWA are displayed 
in Fig. 5 for z — 10 and z = 20, respectively. While it 
is promising for all the three experiments to detect the 
rcionization signals in the angular power spectra of 21cm 
fluctuations over a wide range of redshifts beyond z — 6 
and angular scales from I ~ 10 2 to ~ 10 4 , after an in- 
tegration time of ~ 1 year, a significant detection of the 
BAO wiggles on the 21cm angular power spectrum turns 
to be still difficult especially at higher redshifts (or lower 
frequencies) due to both the weak signals of BAO them- 
selves and the limitation of angular resolutions of cur- 
rent radio telescopes. The strategy is that the maximum 
variations of BAO (e.g. the power difference between 
acoustic peaks and adjacent troughs) should exceed the 
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Fig. 5. — The expected angular power spectra and mea- 
surement errors with 21CMA(top), LOFAR(core) (middle) and 
MWA (bottom) for a bandwidth of Av = 0.1 MHz and an inte- 
gration time of 1 year. 

errors at the corresponding modes. To optimize the de- 
tection, one should choose to work at higher frequen- 
cies near 200MHz, though the total 21cm signals may 
become weaker because the reionization process almost 
completed by z = 6. To be specific, the relatively smaller 
sky coverage and higher system temperature of 21CMA 
prevent a significant detection of the BAO wiggles on 
the angular power spectra of 21cm fluctuations unless a 
longer integration time of ~ 10 years is allowed. In con- 
trast, both LOFAR and MWA may be able to capture 
the BAO signals within ~ 1 year. In particular, MWA 
can even trace all the BAO wiggles out to z w 30 and on 
very small angular scales of a few arcminutes because of 
its numerous independent baselines. We anticipate that 
a similar result can be reached when all the 77 stations 
in the planned LOFAR start to operate. 

5. CONCLUSIONS 

The BAO should be imprinted onto the 21cm emission 
background from the epoch of reionization through the 
underlying density perturbations. Detection of the sig- 
nals will provide valuable information about the forma- 
tion and evolution of cosmic structures at higher redshifts 
beyond z w 6. It also furnishes a standard ruler to the 
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probe of topology and geometry of the universe includ- 
ing dark energy and its equation of state. In particular, 
many of the BAO modes were not erased by the forma- 
tion of large-scale structures by z = 6, and we may be 
able to see the BAO wiggles at smaller scales. This will 
complement our knowledge of the BAO at intermediate 
redshifts in addition to the detections of BAO signatures 
in CMB at z « 1000 and in large-scale distribution of 
galaxies at z « 0. A combination of these BAO mea- 
surements at different redshifts will allow us to set more 
robust constraints on the determinations of cosmological 
parameters. 

We have used an analytic approach based on halo 
model for distribution and evolution of dark matter, in 
which baryons trace essentially dark matter but the stan- 
dard CDM power spectrum is modified by the presence 
of baryons. We have then calculated the 21cm emission 
power spectrum from neutral hydrogen surrounding the 
ionized bubble of each halo, following a simple model 
of ionization history. Our results show that BAO are 
indeed presented at the power spectra of the redshiftcd 
21cm emission from the epoch of reionization, and are al- 
most unaffected by the presence of nonlinear structures 
beyond z > 6. This indicates that one should be able 
to see many of the BAO modes at the 21cm power spec- 
trum. 

We have worked with the angular power spectra of 
the 21cm emission from the epoch of reionization for a 
fixed frequency. The BAO signatures are clearly seen 
at i 500-3000 through the entire history of reioniza- 
tion. However, detections of these wiggles with existing 
and planned radio interferometric arrays such as 21CMA, 
LOFAR and MWA does pose a technique challenge. The 
primary difficulty, apart from the extremely faint signals 
of 21cm emission themselves from the epoch of reioniza- 
tion and the strong foreground at low frequency, arises 
from the high system noise and the limitation of angular 
resolutions. Yet, within an integration time of about 1 
year, it seems that both LOFAR and MWA are capable 
of capturing the BAO signatures at I ~ 1000, provided 
that the foreground contamination can be successfully 
removed. 
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